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Abstract 

We consider integrable Hamiltonian systems in R 2n with integrals 
t^j- | of motion F = (Fi, . . . , F n ) in involution. Nondegenerate singular- 

■ ities of corank one are critical points of F where rank dF = n — 1 

and which have definite linear stability. The set of corank-one non- 
degenerate singularities is a codimension-two symplectic submanifold 
• invariant under the flow. We show that the Maslov index of a closed 

. curve is a sum of contributions ±2 from the nondegenerate singulari- 

ties it encloses, the sign depending on the local orientation and stabil- 
ity at the singularities. For one- freedom systems this corresponds to 
the well-known formula for the Poincare index of a closed curve as the 
oriented difference between the number of elliptic and hyperbolic fixed 
points enclosed. We also obtain a formula for the Liapunov exponent 
of invariant (n — l)-dimensional tori in the nondegenerate singular 
set. Examples include rotationally symmetric n-freedom Hamiltoni- 
ans, while an application to the periodic Toda chain is described in a 
companion paper [TU] , 

1 Introduction 

Maslov indices are integers associated with curves on Lagrangian submani- 
folds of cartesian phase space M 2n = {(q, p)} which count caustics, ie points 
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along the curve where the projection of the Lagrangian submanifold to the 
q-plane becomes singular. Caustics are counted with a sign, which is de- 
termined by the local orientation of the curve near the caustic. For generic 
curves, the Maslov index is invariant under deformations which leave the 
endpoints fixed. In semiclassical approximations to the Schrodinger equation 
(and, more generally, in the short-wave asymptotics of linear wave equations), 
Maslov indices appear as n/2 phase shifts along classical orbits associated 
with focusing of families of orbits. 

In classically integrable systems, Maslov indices appear in the EBK quan- 
tisation conditions fSl ITKj . according to which energy levels of the corre- 
sponding quantum system are given asymptotically by quantising the actions 
variables Ij according to 

h = ( n j + /V 4 )^ j = l,...,n. (1.1) 

fij is the Maslov index of the closed curve (on a Lagrangian n-torus) traced 
out by the conjugate angle variable 9j, and is even in this case. For closed 
curves, Arnold pQ has given a canonically invariant description of the Maslov 
index as a one-dimensional characteristic class, along with an explicit formula 
in terms of a winding number in the space of Lagrangian planes. 

In the last 20 years there has been much interest in the topology of com- 
pletely integrable finite-dimensional Hamiltonian systems. The generic local 
structure and dynamics is given by the Liouville- Arnold theorem [2], ac- 
cording to which neighbourhoods of phase space are foliated into invariant 
Lagrangian submanifolds diffeomorphic to M. n ~ k x T k , where T k is the fc-torus 
and the dynamics is linearised by action-angle coordinates. This local be- 
haviour breaks down at critical points of the energy-momentum map, where 
the dimension of the critical components of the invariant sets drops. A Morse 
theory for integrable Hamiltonian systems, wherein the global topology is de- 
scribed in terms of these critical sets, has been developed by Fomenko jH], 
Eliasson (Zj, Vey j2H| and Tien Zung j2S], among others. The effect of sin- 
gularities on quantum wavefunctions are analysed in Obstructions to 
global action-angle variables were described by Duistermaat jS], and the as- 
sociated phenomenon of monodromy has been studied in the classical and 
quantum contexts 151 12TJ 120] . 

This paper is also concerned with the topology of completely integrable 
systems; nonzero Maslov indices are manifestations of the global topology. 
We obtain a formula for the Maslov index of a closed curve in terms of 
the nondegenerate singularities of codimension two (in M 2n ) enclosed by the 
curve. Our result is analogous to the formula for the Poincare index of closed 
orbits in planar systems in terms of the nondegenerate fixed points enclosed 
by the orbit, and, indeed, for one-freedom systems coincides with it. 
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The relationship between the Maslov index and singularities in the La- 
grangian foliation has been studied from a more general point of view by 
Suzuki j22j. Suzuki considers Lagrangian subbundles of a general symplectic 
vector bundle, (ie, not just the tangent bundle of a cotangent bundle), and 
considers higher Maslov classes of dimension 4k — 3, though does not consider 
the case of integrable Hamiltonian systems. Trofimov [21], who also consid- 
ers higher Maslov classes, shows that the Maslov indices of Liouville tori are 
constant in a connected component of the regular values of the integrals of 
motion. 

Let us fix some notations and conventions. We consider 2n-dimensional 
cartesian phase space M 2n , with canonical coordinates z = (q, p), where 
q, p G M. n . The symplectic inner product of z 1; z 2 G M. 2n is given by 

[zi, z 2 ] = qi ■ p 2 - Pi • q 2 = z 2 ■ J^ 1 • zi, (1.2) 

where 



I 
-I 



1.3) 



and I is the nxn identity matrix. Given a linear subspace A C M 2n , denote its 
skew-orthogonal complement by A z ; this is the subspace of vectors u G M 2n 
for which [u, v] = for all v G A. A is said to be isotropic if A C A z ; in 
this case dim A < n. A is said to be a Lagrangian plane if A z = A; in this 
case dim A = n. The tangent space T z R 2n is naturally identified with IR 2n , in 
which case the symplectic inner product (jl.2|) coincides with the canonical 
symplectic form on T*R n = R 2n . 

Given smooth functions F and G, their Poisson bracket is given by 

{F, G} = dF ■ J ■ dG. (1.4) 

An integrable system on M 2n is described by a Hamiltonian 

H = h(F) (1.5) 

expressed as a smooth function of n smooth, functionally independent func- 
tions F(z) = (Fi(z), . . . F n (z)) in involution, so that 



{F a ,Fp} = 0, l<a,P<n. (1.6) 



Let 



£ a = J-dF a . (1.7) 

£ Q is the Hamiltonian vector field generated by F a . The involution condition 
(jl.fij) is equivalent to 

[Z at tf,] = 0, l<a,(3<n. (1.8) 
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Let 

A = span {£!,...,£„}. (1.9) 

The condition (jl.8j) implies that A(z) is isotropic. If £i(z), . . . , £ n (z) are 
linearly independent, then A(z) is a Lagrangian plane. 

The regular set of the integrable system, denoted R, is the set of regu- 
lar points of F, ie the set on which dFi, . . . , dF n are linearly independent; 
equivalently, R is the set on which A is Lagrangian. R is open in M. 2n . The 
complement of R, the singular set, denoted S, is the set of critical points of 
F. S is the disjoint union of sets k = 1, . . . , n, on which dF has corank 
k, 

S fc = {z | corank dF = k}. (1.10) 

Our interest here is in the set Si, where there is precisely one linear 
relation amongst the dF^s. Given y e E 1; let 



5> Q (y)cZF Q (y) = 0, (l-li; 



Q = l 



where c(y) ^ 0. The relation (jl.ll|) determines c(y) up to a nonzero scalar 
factor. We note that c(y) need not be continuous on Si. ( (|1.11|) defines line 
bundles over submanifolds contained in Si, which need not be trivial; c(y) 
defines local sections of these bundles.) Given y G Si, let 

n 

K(y) = $> a (y)JF Q "(y), (1.12) 

a=l 

where 

d 2 F 
ozoz 

denotes the Hessian of F a . Like c(y), K(y) is determined up to a nonzero 
scalar factor, and need not be continuous on Si. K(y) is an infinitesimal 
symplectic matrix, ie 

[K-u,v] + [u,K- v] = (1.14) 

for all u, v G R 2n . Let 

r=iTrK 2 . (1.15) 

Define the corank-one nondegenerate singular set, denoted A, to be the subset 
of Si on which r ^ 0; 

A = {x e Si |r(x) ^ 0}. (1.16) 
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Note that, as r is determined up to a positive scalar factor, sgnr is well 
defined on A. r(x) is positive (resp. negative) according to whether x is 
an elliptic (resp. hyperbolic) fixed point of c • F (modulo linearly neutral 
directions associated with the action of the remaining components of F) . 

In Section we establish that A is an invariant codimension-two sym- 
plectic submanifold which is invariant under the flow of H. While this follows 
from general results, eg [21], we give an independent argument which serves 
the developments in subsequent sections. In Section we show that the 
Maslov index of a closed curve C in the regular set is typically given by a 
sum of contributions ±2 from the nondegenerate singular points it encloses, 
the sign determined by the local orientations and signatures sgnr at the 
singularities. (Typical cases include those for which £ — A is contained in a 
submanifold of codimension at most three.) In SectionHJ we show that sgnr 
is related to the linear stability of A; given x e A with compact orbit under 
the integrable flows, we derive a formula for the Liapunov exponent, which 
vanishes or not according to whether r(x) is positive or negative. Thus, the 
signs in the Maslov index formula may be regarded as a product of signs 
determining the orientation and stability of the singularities. 

Some simple examples are discussed in Section EJ For one-freedom sys- 
tems, our result coincides with the well-known formula for the Poincare index 
of a closed curve C in a planar vector field. Examples where the formula does 
not apply, which can arise in bifurcations, are also discussed. In Section |B] 
we consider rotationally invariant systems in W 1 . In a companion paper, we 
consider a non-separable example, the periodic Toda chain [TO] . 

In what follows, we use the notation A = B to indicate that A and B 
differ by a nonzero scalar factor (ie, they are projectively equivalent). A and 
B may be nonzero scalars, vectors, or matrices. 

2 The nondegenerate singular set 

We show that the nondegenerate singular set A is a codimension-two sym- 
plectic submanifold invariant under the integrable flow. In outline, the ar- 
gument is as follows: We introduce a complex-matrix- valued function, M(z), 
whose determinant vanishes precisely on X. It is shown that A is an open 
subset of the set of regular points of det M = 0, and a calculation establishes 
that the symplectic form is nondegenerate on A. For brevity, we will denote 
det M by |M|; thus A is characterised by |M| = 0. 

It turns out that the Maslov index of a closed curve is twice the winding 
number of arg |M| along the curve. This fact, along with the expression ()2.5)1 
for d|M| derived below, is the basis for the Maslov index formula derived in 
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Section EJ 

Let M(z) be given by 



The a th column of M has real and imaginary parts equal to the q- and (— p)- 
components of £ Q respectively. Then the set of corank-/c singularities E*., as 
given by (jl.lUj) . may be characterised as follows: 

Proposition 2.1. z 6 Sj. -<=>- corankM(z) = k. 

Proof. If X)a=i a ctdF Q (z) = for some nonzero a £ IR n , then it is straightfor- 
ward to show that M(z) - a = 0. Therefore, dim span {dF i(z), . . . , dF n (z)} > 
rankM(z). On the other hand, suppose a £ C n is a nonzero right nullvec- 
tor of M(z). Then a is also a right nullvector of (M^M)(z). The involution 
condition (jl.6|) implies that (MtM)(z) is real. Therefore, without loss of 
generality, we may assume that a is real. For a real, it is straightforward 
to show that M(z) • a = implies that Y^=i a adF a (z) = 0. Therefore, 
rankM(z) > dim span {dF x {z), . . . ,dF n (z)}. □ 

For y £ Ei, the preceding implies that c(y), as given by (jl.llj) . spans the 
right nullspace of M(y). Let b(y) £ C n be the corresponding (complex) left 
nullvector, ie 

b(y) • M(y) = 0. (2.2) 

Like c(y), b(y) is determined up to a (complex) nonzero scalar factor, and 
need not be continuous on Si. 

For y £ Ei, let (3(y) £ M 2n be given by 

/3(y) = (Reb(y),Imb(y)). (2.3) 

Let 

r/(y) = (K-/3)(y), 

0(y) = (KJ-/3)(y), (2.4) 

where K(y) is given by ()1.12jl . Then we have the following formula for the 
derivative of the determinant of M: 

Proposition 2.2. d|M| vanishes on E^ for k > 1. For y £ Ei 7 

d|M|(y)-J-(77(y) + i0(y)). (2.5) 
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Proof. We have the general formula for the derivative of a determinant, 

d|M| = Tr (m T dM) , (2.6) 

where m(z) is the cofactor matrix of M(z). From Proposition 12.11 if k > 1, 
then m(z) vanishes on Therefore, d|M| = on for k > 1. 

Let yo G Si. Choose fixed matrices So and To with unit determinant 
whose first columns are c(y ) and b(y ), respectively. Let N(z) = TqM(z)So- 
Then |N|(z) = |M|(z), and 

d|M| = Tr(n T dN), (2.7) 

where n(z) is the cofactor matrix of N(z). By construction, the first row 
and first column of N(yo) vanish, so that n(yo) has a single nonzero element, 
namely rin(yo). It follows that 

d|M|(y„) = n u (y )diVii(yo) = dN n (y ) = (b • dM ■ c) (y ). (2.8) 

Some straightforward manipulation shows that the preceding is equivalent to 
(J23J) above. □ 

From Proposition 12. 2\ it follows that the regular component of the level 
set |M| = 0, ie, the subset of the level set on which d\M\ has maximal 
(real) rank equal to two, is contained in Si, and is the set where r\ and 6 
are linearly independent. Linear independence is implied by the stronger 
condition [77, 9] ^ 0, which, in turn, is equivalent to the following condition 
on K: 

Proposition 2.3. Let y G Si ; and let kerK 2 (y) and imK 2 (y) denote the 
kernel and image o/K 2 (y) respectively. If r(y) 7^ ; then dimimK 2 (y) = 2, 
dimkerK 2 (y) = 2n — 2, and 

ker K 2 (y) © im K 2 (y) = R 2n (2.9) 

is a decomposition of M 2n into symplectic skew-orthogonal subspaces. More- 
over, 

imK 2 (y) =span{r7(y),0(y)}, (2.10) 

and 

sgn[r/(y),0(y)] = sgnr(y). (2.11) 

Proposition 12. 31 can be deduced from Williamson's theorem (see, eg, [2]). We 
give an explicit argument below. One can also show that, if r(y) = 0, then 
K 2 (y) = 0and [r/(y), 0(y)] = 0. 
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Proof. Let E denote the generalised nullspace of K(y) (the nullspace of 
powers of K(y)) and let r be a positive integer such that K r (y) • E — 0. Let 
E* = im K r (y). Then 

E ®E, = R 2n . (2.12) 

The fact that K(y) is infinitesimal symplectic implies that Eq and E* are 
skew-orthogonal. For if u G Eq and v = K r (y) ■ w G E*, then 

[u, v] = [u,K r (y) • w] = (-l) r [K r (y) ■ u, w] = 0. (2.13) 

As Eq and E* are complementary and skew-orthogonal, the restriction of 
the symplectic inner product to either is non-degenerate, so that both are 
symplectic subspaces of M 2n . 

Next, we show that K 2 (y) • Eq = 0, so that we may take r = 2 above. We 
have that 

K(y)-^(y) = 0; (2.14) 

this follows from differentiating {F a , Fp} = 0, multiplying the resulting equa- 
tion by c Q (y) and summing over a. Thus K(y) ■ A(y) = 0, so A(y) is an 
isotropic subspace of E . Since y G Si, dim A(y) = n — 1. As E is symplec- 
tic, it follows that dimE is either 2n — 2 or 2n. But the latter would imply 
that TrK 2 (y) = 2r(y) = 0, contrary to assumption. Thus dimi^o = 2n — 2 
and dimK* = 2. 

The fact that dimi^o = 2n — 2 implies that A(y) is a maximal isotropic 
subspace of E ; ie A z (y) n E = A(y). It follows that K(y) ■ E C A(y). For 
if u G Eq and v G A(y), then [K(y) • u, v] = — [u, K(y) • v] = 0, so that 
that K(y) • Eq is skew-orthogonal to A(y). Since A(y) is a maximal isotropic 
subspace of Eq, it follows that K(y) • Eq C A(y). Thus, K 2 (y) • Eq = 0, so 
that 

£ = kerK 2 (y), E t =imK 2 (y). (2.15) 

Let K* denote the restriction of K(y) to E*. From (|2.12jl . TrK* = 
TrK(y) = (K is necessarily traceless) and TrK 2 = TrK 2 (y) = 2r(y) ^ 0. 
It follows that the characteristic polynomial of K* is given by 

K^-r(y) = 0. (2.16) 

Let 

(3 = (3 + f3„ J./3 = 7o + 7 , (2.17) 

denote the decompositions of (3 and J • (3 into their respective components in 
Eq and E*. From ()2.2j) and (|2.Hj) . one can verify that /3(y) and J • /3(y) are 
skew-orthogonal to A(y). Since f3^ and 7^ are necessarily skew-orthogonal 
to A(y), it follows that (3 and 7 are as well, so that 

/3 , 7o e^onA z (y) = A(y). (2.18) 



Then 

(3 2 = [J ■ (3,(3} = [ 7 „/3J + [ 7o ,/3 ] = fr*,P*]- (2-19) 

Now let r/(y) = K(y) ■ /3(y), 0(y) = K(y) • 7 (y), as in flQJ. From d2H2J 
and pTRjl . 

I7(y) = K, ■ /3„ 0(y) = K* • % . (2.20) 
Therefore, from (|2~THjl . (j2~lHJ) and JZZfl) , 

fo(y),0(y)] = [K* • (3*, K* • 7 J = -[K 2 • j8„ 7 J = -r(y)[/3„ 7 „] = r(y)/3 2 , 

(2.21) 

from which ()2.11|) follows. This in turn implies that 7](y) and 0(y) are lin- 
early independent, so that, from (|2.2()jl . imK 2 (y) = E# — span{?7(y), 0(y)}. 

□ 

Let $t(z) denote the flow of the Hamiltonian H, and let S(z, t) denote 
the linearised flow, ie the 2n-dimensional symplectic matrix given by 

Sz,t ) = — * z . 2.22 
oz 

S(z, i) satisfies the differential equation 

S(z,t) = JH /, ($t(z))S(z,t). (2.23) 

The following shows that Si is invariant under the flow, while the quantities 
c and K are invariant up to a nonzero scalar factor. 

Proposition 2.4. IfyE Si, i/ien y t = 3>t(y) G Si, and 

c(y t ) = c(y) (2.24) 
K(y t ) = S(y,t)K(y)S- 1 (y,t). (2.25) 

Proof. Integrability implies that the vector fields £ Q (z) are invariant under 
the flow, so that 

=S(M) •*<■(*)■ (2.26) 

Since S(z,t) is invertible, dimA(z) is invariant under the flow, which implies 
that the regular set R and the components S^ of the singular set are sep- 
arately invariant. Let y G Si. (|2.26J) implies that Yla=i c a(y)$ a {yt) = 0. 
As y t G Si, there is just one linear relation amongst the £ a (yt), so ()2.24|) 
follows. ()2.25j) follows from differentiating (|2.2fij) to get 

JF/(^(z)))S(z, t) = ^(z, t) ■ U*) + S(z, t) JF a "(z). (2.27) 
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Letting z = y in the above, multiplying by c a (y) and summing over a, we 
get 

n 

5> Q (y)JF a "( yi ) = S(y,t)K(y)S" 1 (y,t). (2.28) 

Ol=l 

From ()1.12|) and (|2.24|) . the left-hand side of the preceding is, up to nonzero 
scalar factor, K(y t ). (J2.2f>j) follows. □ 

Definition 2.1. The corank-one nondegenerate singular set, denoted A, is 
the subset o/Si on which 

Theorem 2.1. A is a codimension-two symplectic submanifold which is in- 
variant under the flow. Given x G A, T X A = kerK 2 (x) and (T X A) Z = 
imK 2 (x). 

Proof. The regular component of the level set Si = {z| |M|(z) = 0}, ie the 
subset of Si on which rf|M| has maximal rank, is a (2n — 2)-dimensional 
submanifold. From Propositions 12.21 and 12. 3| A is contained in the regular 
component of Si. Since A is open in Si, it follows that A is itself a (2n — 2)- 
dimensional submanifold. From Proposition l2.2l the tangent space T X A is the 
skew-orthogonal complement of span {rf(x), 0(x)}. From Proposition 12.31 it 
follows that T X A = ker K 2 (x) and (T X A) Z = imK 2 (x). From Proposition l2.4[ 
Si is invariant under the flow, and for y G Si, r(y 4 ) = r(y). It follows that 
A is invariant. □ 



3 Singularity formula for the Maslov index 

If z belongs to the regular component R of an integrable system, then A(z) 
is a Lagrangian plane. Arnold £Q showed that A(n), the space of Lagrangian 
planes, has fundamental group 7i"i(A(n)) = Z, with continuous closed curves 
characterised by an integer winding number. 

Let C denote a continuous, oriented closed curve in R, parameterised by 
z (s), < s < 1, with z(l) = z(0). Then A(C) describes a continuous, ori- 
ented closed curve in A(n), parameterised by A(z(s)). We define the Maslov 
index of C, denoted /i(C), to be the winding number of A(C) in A(n), ie 

fi{C) = wnA(C). (3.1) 

Arnold |]Q showed that, under certain genericity conditions, for curves C on 
a Lagrangian manifold (eg, an invariant torus of an integrable system), the 
definition (|3.1|) coincides with the signed count of caustics along C (for curves 



10 



on invariant tori, this is the Maslov index which enters into the semiclassical 
quantisation conditions (jl.ljO . 

We shall make use of the following explicit formula for the Maslov index 

[IB!: 

/i(C) = i(arg|M|(z(l))-arg|M|(z(0))), (3.2) 

where M is given by (|2.1|) . and arg |M|(z(s)) is taken to be continuous in s 
(this is possible because, from Proposition 12. 1[ |M| does not vanish in R). 
Thus, fi{C) is twice the winding number of the phase of |M| evaluated along 
z(s). (The formula ()3.2j) may be obtained by averaging the count of caustics 
with respect to a one-parameter family of projections. A related formula is 
given in [14] . and an application to resonant tori is discussed in |19j.) 

If C is contractible in R, then A(C) is contractible in A(n), so that n(C) 
must vanish. Therefore, if C has a nonzero Maslov index, it must enclose 
points in the singular set. Typically, the Maslov index can be determined by 
the singularities enclosed. This is most easily demonstrated for a small curve 
about x G A in the plane spanned by (T X A) Z . Let C° denote the family of 
curves 

z°( s ) = x + e(cos 2ns tj(x) + sin 2ns 0(x)). (3.3) 

For sufficiently small e ^ 0, C° is contained in R, and the only singular point 
it encloses is x (that is, C° can be contracted to x without passing through 
any other singular points). From Propositions 12. II and 12 ,2| 

|M|(z°(s)) = |M|(x) + ed|M|(x) • (cos 2tts 77(x) + sin 2vrs 0(x)) + 0(e 2 ) 
= const x iefa(x), 0(x)]e 2ms + 0(e 2 ), (3.4) 

where const denotes a real, nonzero constant independent of e, and [^(x), 0(x)] ^ 
by Proposition 12.31 Thus, 

arg|M|(z°(s)) - arg |M|(z°(0)) = 2ns + 0(e 2 ). (3.5) 

It follows from ()3.2j) that for sufficiently small e, 

MC e °) = 2. (3.6) 

More generally, consider the family of curves C e given by 

z e (s) = x + e(cos27rs u + sin 2ns v), (3.7) 

where u, v G T x IR 2n . Then 

|M|(z £ (s)) = |M|(x) + ed|M|(x) • (cos2ttsu* + sin27rsv*) + 0(e 2 ), (3.8) 
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where u*, v* denote the projections of u, v to (T X A) Z . If u* and v* are lin- 
early independent, then, for sufficiently small e, /x(C e ) = +/i(C°) or — /u(C°) 
according to whether C e and C° are similarly or oppositely oriented, ie ac- 
cording to whether sgn [u*, v*] = ±[?7(x), 0(x)]. Since sgn [t/(x), 0(x)] = 
sgnr(x) (cf Proposition 12. 3j) . we get that 

fi(C e ) = 2 sgn [u*, v„] sgnr(x). (3.9) 

We can extend ()3.9j) to larger curves in R as follows. Let D 2 = {(x, y) \x 2 + 
y 2 < 1} denote the unit two-disk endowed with its standard orientation. 

Definition 3.1. A nondegenerately transverse disk is a continuous map 
S : D 2 — > IR 2ri smooth on the interior of D 2 such that i) the image of S 
is contained in the union i?UA of the regular set and the corank-one nonde- 
generate singular set, ii) S~ 1 (A) ; the preimage of the singular set, is a finite 
set of points tj = (xj,yj) in the interior of D 2 , and Hi) n*odS(ej) : T ej D 2 —> 
(T x . A) z ; where Xj = S(ej) and 7r* is the projection onto (T x . A) z with respect 
to the symplectic decomposition M. 2n = T K . A © (T x A) z , is nonsingular. 

That is, S is a nondegenerately transverse disk if its image intersects S 
transversally at a finite set of points in A. We define aj, the local orienta- 
tion of Sj at Xj, to be +1 (resp. —1) if 7r* o dS(ej) is orientation-preserving 
(resp. orientation- reversing), with the orientation on (T X .A) Z determined by 
the restriction of the symplectic form. Explicitly (cf (|2.16j) ). 

(jj = sgn [Uj>, Vj-J = sgnr(xj) sgn [K 2 (xj) • u,-, vj , (3.10) 

where 

dS, , dS . . 

and Uj*, Vj* denote the projections of u, v in (T x . A) z . 

Theorem 3.1. Lei C G R be the oriented boundary of a nondegenerately 
transverse disk S, with C parameterised by z(s) = S(cos27rs,sin27rs). T/ien 

/i(C , )=2^a J sgnr(x J ). (3.12) 

5=1 

where the sum is taken over preimages ej G S _1 (A), and x^ = S(ej). 

If the set E — A is contained in a submanifold of codimension three or 
more, then it can be shown that any closed curve in R can be realised as the 
boundary of some nondegenerately transverse disk (the argument is omitted). 
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In this case, (|3.12|) applies to all C C R. If E — A is of codimension two, then 
()3.12j) may not apply. Such cases can arise in connection with bifurcations, 
as discussed in Sectional 

Theorem 13.11 may be summarised as follows: if C encloses only isolated 
nondegenerate singular points, then its Maslov index is a sum of contributions 
±2 from each of the singularities Xj enclosed, the signs depending on the 
local orientations <jj and sgnr(xj) at the singularity. In Section 01 it is 
shown that if the orbit of Xj under each of the integrable flows is compact, 
then sgnr(xj) = —1 if these orbits are linearly stable, and sgnr(xj) = +1 if 
they are linearly unstable. Thus, sgnr(xj) is determined by the stability of 
the singularity. We note that if S s is a continuous family of nondegenerately 
transverse disks with fixed boundary C, both the number of singularities and 
their local orientations and stabilities may vary with s. However, the signed 
sum of the products of their orientations and stabilities remains invariant. 

Proof of Theorem \3.1\ Choose r < 1 so that all of the singular preimages 
Cj G D 2 are contained in the disk x 2 + y 2 < r 2 . By continuity, the Maslov 
index of C is equal to the Maslov index of the image of this circle under S. 
Since S is smooth on the interior of D 2 , it follows from ()3.2j) that 

n(C) = -<f darg|M|oS. (3.13) 

By Stokes' theorem, the integration contour can be replaced by a sum of N 
positively oriented circles centred at each of the e/s of radius e, with e taken 
to be small enough so that the circles do not overlap. Let Cj C R denote the 
images of these circles under Sj. Then 

N 

/*(<?)= X>(Q)- (3.14) 

i=i 

Cj is parameterised by Zj(s) given by 

Zj(s) = Xj + e(cos 2irsUj + sin27rsVj) + 0(e 2 ), (3.15) 

where Uj and Vj are given by (|3.11|) . For e sufficiently small, the 0(e 2 ) 
terms can be dropped without changing the Maslov index of Cj. Since Sj 
is nondegenerately transverse, the projections of Uj and Vj to (T X A) Z are 
linearly independent. From ()3.9j) and (|3.1()j) . 

fi(Cj) = 2a i sgnr(x). (3.16) 

Substituting into (|3.14|) . we get the formula (|3.12j) . □ 
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4 Transverse linear stability of the nondegen- 
erate singular set 

Let x be a nondegenerate corank-one critical point which is fixed by the flow 
of c-F. Clearly, the linear stability of x under c-F is determined by the spec- 
trum of K(x); in particular, x is elliptic or hyperbolic according to whether 
t(x) is negative or positive. Here we consider the linear stability of x under 
the Hamiltonian H = h(F), which need not leave x fixed. We assume that 
the F-orbit of x is compact. One way to address this question would be 
through a systematic normal form description of the dynamics in a neigh- 
bourhood of the F-orbit of x (see j^E] and, for more detailed formulations, 
jH] for two degrees of freedom and ^Hj for n degrees of freedom). Instead, 
we here obtain directly an explicit formula for the Liapunov exponent of the 
if -orbit of x. 

Let ( z ) denote the Hamiltonian flow generated by F a (z), a = 1, . . . , n. 
Integr ability implies that these flows commute. Given x £ A, let 



where s = (si, . . . , s n ) £ W 1 . Then x i— > x s defines an IR n -action on A. 
Let T x denote the orbit of x under this action. The identity component 
of the isotropy subgroup of x is the ray s = sc in M. n , where c • dF(x) = 
0. Quotienting out by this subgroup gives an M n_1 -action on T x . If 7^ is 
compact, the Liouville- Arnold theorem implies that it is topologically an 
(n — l)-dimensional torus, and with Theorem 12 . II that the Hamiltonian flow 
generated by H describes a (2n — 2)-dimensional integrable system in a 
A-neighbourhood of T x . In what follows we assume that T x is compact, and 
let (■)%, denote the average over 7^ with respect to the normalised invariant 
measure. The linear stability of x is determined by the (maximal) Liapunov 
exponent ft#(x), given by 



where S, the linearised flow, is given by (|2.22|) . and 1 1 ■ 1 1 denotes the Euclidean 
norm on M 2n . It is a standard result that the Liapunov exponent exists for 
almost all initial conditions and is independent of the choice of metric (see, 



x s = <(^ 2 (- ••*:>))•••), 



(4.1) 




(4.2) 



eg, [12]). 



For 



x s £ %t given by (|4.1j) . let 



n 




(4.3) 



a=l 
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Clearly, K(x) = K(x), and from Proposition 12.41 f which holds in particular 
for H = F a , a — 1, . . . , n), it follows that 

K(x s )-K(x s ). (4.4) 

Unlike K, K is necessarily smooth on T x . Arguing as in (|2.27[) . we have that 

K(x s ) = U(x, s)K(x)U" 1 (x, s), (4.5) 

where 

U(x,s) = ^. (4.6) 

It follows that K(x s ) has a pair of nonzero eigenvalues ±r 1 / 2 (x). If r(x) > 0, 
Proposition 12.31 implies that 

Q(x s ) = (K(x s ) + r 1 /2( x ))K 2 (x s ) (4.7) 

is, up to normalisation, a projector onto the one- dimensional r 1 / 2 (x)-eigenspace 
of K(x s ). In general, Q(x s ) is not symmetric. The normalised symmetric 
projector onto the one- dimensional r 1//2 (x)-eigenspace is given by 

p <-> = (t^q) W- < 4 ' 8 > 

Theorem 4.1. Let x G A, and suppose that T x is compact. If r(x) < ; 
then kh(x) = 0. 7/r(x) > ; then 



j2(^r(Fi,...,F n )\x) Ka (x), (4.9) 

a=l a 



where K a (x), the Liapunov exponent for H = F a , is given by 

« a (x) = |<Tr(PJF a ")> r J, (4.10) 

and P is given by fl^.£| ). The Liapunov exponents « a (x) do not all vanish; in 
particular, 

n 

^c q (x) Kq (x) = t 1 / 2 (x). (4.11) 

a=l 

We note that the transpose Q T in ()4.8j) is defined with respect to the Eu- 
clidean inner product. For a general non-Euclidean metric one would a obtain 
a more general expression for K a ; however, its value would be unchanged. 
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Proof. As noted above, in a neighbourhood of 7^, the restriction of the flow of 
H to A, regarded as a symplectic manifold of dimension 2n — 2, is integrable. 
As Liapunov exponents for compact integrable flows vanish, the right-hand 
side of (14. 2j) vanishes for x tangent to A, so we may restrict x to the two- 
dimensional transverse plane £*(x) defined in (|2.15|) . As nonzero Liapunov 
exponents for Hamiltonian systems occur in signed pairs, for any nonzero 
X £ the limit on the right-hand side of ()4.2|) either vanishes, or else is 

equal to ±%(x). Therefore, if one takes the absolute value of the expression 
on the rhs of (J4.2)) . the supremum over % is no longer necessary. 

We consider first the case r(x) = — uj 2 < 0. For x s £ T x , K(x s ) has a 
pair of imaginary eigenvalues ±za>. Let C(x s ), C*( x s) denote corresponding 
conjugate eigenvectors. The real and imaginary parts of C( x s) span i?»(x s ). 
Therefore, [C( x s) ; C*( x s)] cannot vanish. The normalisation condition 



determines C( x s) up to a complex phase factor (we do not assume that this 
phase factor can be chosen to make C( x s) continuous on %_). 

Let x< = $t(x). ()4.5j) implies that S(x, t) ■ £(x) is proportional to C( x 0- 
Since the symplectic inner product is preserved under the linearised flow, the 
normalisation condition ()4.12|) implies that S(x, t) ■ C( x ) differs from C( x t) 
by a phase factor, so that 



where | |C| | 2 = 1 1 Re C| | 2 + I l im CI | 2 - Since | |C( x s)| I is bounded on 7*, | |S(x, t) ■ 
C( x )|| is bounded in t, so that /^(x) = 0. 

Next we consider the case r(x) > 0. For convenience, we first assume 
that H = F a . Take x £ -E*( x ) to be an eigenvector of K(x) with eigenvalue 
r 1 / 2 (x). We may write (|4.2j) as 



[C(x s ),C(x s )] = z 



(4.12) 



S(x,t)-C(x)|| = ||C(x,)||, 



(4.13) 




(4.14) 



where xif) — S(x, t) ■ X- From 



d_ 

It 



S(x,t) = JF Q "(x t )S(x,t) 



(4.15) 



it follows that 



d 
It 



log ||x(t)|| 



Xit) ■ JF a w (x t ) ■ xit) 

x(t) ■ x(t) 



Tr (P(x t )JF Q "(x t )) 



(4.16) 
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where P is given by ([4.8)1 . The fact that the commutative flow ()4.1|) is 
transitive on 7^ implies that the Liapunov exponent is constant on T^, and 
therefore is equal to its average. (J4.14)) and (|4.16|) . 



K n x 



(«a(Xs)) 



lim - / T Tr (P((x s ) t )JF/((x s ) t )) 



• (4-17) 



If the integrand in ()4 . 1 7|) is averaged over T x , the average over t becomes 
redundant. We obtain 



n a (x) = |(TrPJF Q "> 



r x | • 



(4.18) 



For general H = h(Fi, . . . , F n ), the expression for the Liapunov exponent 
is obtained by replacing JF a " by JH" in the preceding calculation. We have 
that 

X • JH " • * = E ■ JF a " -X+ £ a^(* ' ° fe' X] ■ (419) 



a=l 



a,/3=l 



For % G -E*(x), the second term vanishes (cf Proposition 12. Hj) . Thus, we get 

dh 



E 



Q = l 



- 



X K„ X 



(4.20) 



in accord with (|4.9jl . If we let H(z) = ^a=i c «-^ Q ( z ) an< ^ use ^ ne ^ ac ^ that 
P Y, n a =i c « JF a = PK = r 1 / 2 ?, we get from (j^pj) that 



E 

a=l 



c a lx)re a (x) = r 1/2 (x 



(4.21) 



as in (jUIH) 



□ 



5 Examples 

5.1 One- freedom systems. 

Let H(q,p) be a smooth Hamiltonian on IR 2 . The singular set E consists of 
critical points of H, ie fixed points of the flow, and the nondegenerate singular 
set A consists of isolated critical points where Tr (JH") 2 = —2det(JH") ^ 0. 
These are the hyperbolic (det JH" < 0) and elliptic (det JH" > 0) fixed points 
of H. The space of Lagrangian planes A(l) is just the projective line RP 1 , 
and the Maslov index (jH.lj) of a closed oriented curve in R = M 2 — E is 
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just (—2 times) the Poincare index (see, eg, [TT] ) of the velocity field JdH 
around the curve. If all the fixed points are isolated and nondegenerate (ie, 
£ = A), then the result ()3.12j) is equivalent to the standard expression for 
the Poincare index as the number of elliptic fixed points minus the number 
of hyperbolic fixed points enclosed by a (positively oriented) curve. 

The Hamiltonian H = p 2 /2 — q 3 /3 has a single degenerate fixed point at 
the origin. Formula ()3.12|) does not apply in this case, but it is straightfor- 
ward to show that the Poincare index (and Maslov index) vanishes for every 
closed curve in the plane, whether or not it encloses the origin. The Hamilto- 
nian Hq can be embedded in a one-parameter family H a = p 2 /2 — (q 3 /3 + aq) 
which undergoes a saddle-centre bifurcation at a = 0, and for which ()3.12|) 
applies for a ^ 0. The index of a closed curve C about the origin is constant 
through the bifurcation, vanishing for a — > 0— because there are no fixed 
points, and for a — > 0+ because the contributions from the saddle at —\fa 
and the centre at i/a cancel. A two-dimensional example is given next. 

5.2 Two-freedom bifurcation. 

The integrable system 

Fi = \v\ + |P2??, F 2 = p 2l (5.1) 

for which q 2 is an ignorable coordinate, is a nongeneric example where the 
degenerate singular set X — A is of codimension two. The singular set is 
given by p\ = 0,p 2 <?i = 0. With appropriate choices of c, TrK 2 = — p 2 . 
The nondegenerate singular set A has two disconnected components, namely 
Pi = <?i = 0,p 2 > (stable fixed point in the (gi,pi)-plane) and Pi = qi = 
0,p 2 < (unstable fixed point in the (gi,pi)-plane.) The degenerate singular 
set is the coordinate plane p = 0. The closed curve C, with q fixed and 
p = e(cos#, sin 9), encloses only the degenerate singular set (see FigureH^b)). 
A calculation, for example using (J3.2J) . shows that /i(C) = 2. 

The degeneracy can be lifted by embedding the system in a one-parameter 
family, 

Ft =b? + |j>2<Z?-«Zi, F 2 =p 2 . (5.2) 

The singular set is given by p\ = 0,p 2 qi = e, and TrK 2 = — p 2 as before. The 
system undergoes a transcritical bifurcation at e = 0; for e ^ 0, the entire 
singular set is nondegenerate. The Maslov index of C is independent of e (at 
least for e small). For e < 0, C encloses a line of stable fixed points with 
positive orientation, as determined by ()3.1()j) . For e > 0, C encloses a line of 
unstable fixed points with negative orientation. See Figure ^a) and^c). 
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Figure 1: Cycle in phase space with H = ^pf + \p2cft — eqi- A*(C) = 2. (a) 
For e < 0, C encloses a stable singularity with positive orientation, (b) For 
e = 0, C encloses a degenerate singularity, (c) For e > 0, C encloses an 
unstable singularity with negative orientation. 



6 Rotationally invariant Hamiltonians in W 1 
6.1 Definitions 

Consider 2n-dimensional Euclidean phase space M. 2n with points denoted z = 
(r, p). Let r(j) G M J denote the projection of r to its first j components, 
similarly p^) the projection of p. Let z^) = (ry), P(j)- Let rh\ = ryj • r^), 
and similarly p^. 

The standard action of the rotation group SO(n) on ~R n , r i— > TZ ■ r, lifts 
to the canonical action (r, p) i— > (7?. • r, 1Z ■ p) on IR 2n . The Hamiltonian 
generators are the angular momenta L a p given by 

L a /3 = r a p/3 - rpp a , 1 < a, (3 < n. (6.1) 

Clearly L a p = —Lp a . The Poisson bracket of components of angular mo- 
menta are given by 

{La/3, L^s} = L ai 5ps — L a s5/3 7 — L^ y S a s + Lp$5 ai . (6.2) 

Let 

4) = \ Ei<«^ L lp = ijfl) - (*W • PO)) 2 ' 2 < j < n. (6.3) 

L 2 n ^ is the squared total angular momentum, which we also denote by L 2 . 
More generally, L 2 ^ is the squared total angular momentum of the projection 
of z to W x W . For future reference, we note that 

L 2 (j) = L 2 U _ 1} + z m ■ Qy_i) ■ z {j] , (6.4) 

where zy] = (rj,Pj) and 

Q«-i,= f ^ ' V' P ' ; V (6.5) 

It is easily verified that 

= if a,/3 < j, {L 2 u} ,L 2 k) } = 0. (6.6) 

Consider a rotationally symmetric Hamiltonian H, characterised by 

{H,L al3 } = 0, l<a,P<n. (6.7) 

Then H is integrable. As integrals of the motion, we may take F\ = H, 
Fi = L±2, and Fj = lJh\ for 3 < j < n. In fact, rotationally symmetric 
Hamiltonians in R n are superintegrable ^7j; they possess (at least) 2n — 2 
independent constants of the motion. It would be interesting to incorporate 
superintegrability into our treatment, but here we will treat if as a standard 
integrable system. 
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6.2 Singularities 

Singularities are points z = (r, p) for which 

n 

d dH + c 2 dL 12 + ^ cj dL 2 {j) = (6.8) 
i=3 

is satisfied for nonzero c e M n . We choose a basis for the c's (by taking 
successive sets of components to vanish). For each basis element, we obtain 
necessary and sufficient conditions for ()6.8|) to hold and for the corresponding 
singularities to be corank-one nondegenerate. 

6.2.1 ci 7^ 0. Spherical singularities. 

If ci ^ in (jfj.Hjl . then at the singularity, the Poisson bracket of H with any 
rotationally invariant function / must vanish, since i) {/, H} = df • J -dH, ii) 
dH may be expressed in terms of the dL a ps (cf (|6.8)l ). and iii) {/, L a p} = df ■ 
J ■ dL a/ 3 = 0. There are three functionally independent rotational invariants, 
eg r 2 , p 2 and r • p. The conditions {H, r 2 } = {H, p 2 } = {H, r • p} = at the 
singularity are equivalent to 

r ■ H p = 0, p • H r = 0, r-H r = p-H p . (6.9) 

The conditions ()6.9)1 also imply a singularity. To see this, note that because 
H is rotationally invariant, it can be expressed as a function of the invariants 
r 2 , p 2 and r • p. It follows that the derivatives of H are of the form 

H r = fr + gp, H p = gr + hp, (6.10) 

where /, g and h are functions of the invariants. Then (|6.9J) implies that 

H r = k(p 2 r — (r • p)p), H p = k(r 2 p — (r • p)r), (6-H) 

where k = f/p 2 = —g/{v ■ p) = h/p 2 . Therefore, 

dH = kdL 2 , (6.12) 

which is just (16. 8J) with c\ = 1, c n = —k, and all other Cj's equal to zero. The 
if-orbits through such points reduce to fixed points of the radial motion. For 
this reason, we call these spherical singularities. 

We determine next the condition for spherical singularities (r, p) to be 
corank-one nondegenerate. We may assume that r and p are not both zero 
(as all of the dF^s vanish at the origin). For definiteness, let us assume 
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that r ^ (the treatment for p ^ is similar). In the neighbourhood 
of the singular point, we introduce local canonical coordinates r, p r , L, 9 
and Z, where p r = p • r/r denotes the radial momentum, L = \fl? the 
magnitude of total angular momentum with conjugate variable 9, and Z the 
remaining 2n — 4 canonical coordinates. As H is rotationally invariant, it can 
be expressed locally as a function of r, p r and L, ie H = h(r,p r ,L). Then 
K 2 has only one nonvanishing two-dimensional block, which corresponds to 
the (r,p r )-plane. Then 

lTrK 2 = -(h rr h PrPr -h* pr ). (6.13) 

Thus, spherical singularities are nondegenerate if the corresponding radial 
fixed points are nondegenerate. 

6.2.2 (03, . . . , c re ) ^ 0, c\ — 0. Axial singularities. 

Let m denote the highest index for which c m ^ 0. We take c m = 1, so that 
()6.8p takes the form 

m— 1 

c 2 dL 12 + Cj dL l) + ^(m) = °- ( 6 - 14 ) 
3=3 

The term dL 2 ^ contains the one- forms dr m and dp m , whereas the other 
terms do not. From ()6.4j) . the condition for the coefficients of dr m and dp m 
to vanish is that 

Q(m-i) • Z[m] = 0. (6.15) 

Let us suppose that 

detQ (m _i)^0 (6.16) 

(the case det Q( m _i) = is considered in Section 16.2.31 below). Then ()6.15|) 
implies that 

Z[ m ] = (r m ,Pm) = 0. (6.17) 

It is easily seen that (J6.17j) also implies a singularity. For if Z[ m \ vanishes, it 
follows from ()6.4|) that 

dL? m _ x) - dL\ m) = (6.18) 
for m > 3, and for m = 3 that 

dLl - 2L 12 dL 12 = 0. (6.19) 
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For both (j6.18|) and (j6.19|) . K has a single nonvanishing two-dimensional 
block equal to JQ( m _i). It follows that 

±TrK 2 = -detJQ (m _D = -L 2 m) . (6.20) 

Thus, the nondegeneracy condition is L 2 m ^ 7^ 0, in which case the singularity 
is elliptic. We call singularities with r m = p m = m- axial singularities, as 
the projections of their if -orbits to the coordinate and momentum m-planes 
in fact lie in the respective (m — l)-planes. 

6.2.3 (03, . . . , c n ) 7^ 0, ci = 0. Radial singularities. 
If 

detQ (m _ 1) = L^ 1) = 0, (6.21) 

then (j6.15j) is satisfied by taking Z[ m ] to be a null vector of Q( m -i). It then 
follows from ()6.4|) that L 2 ^ = 0, and hence 

dL 2 m) = 0, (6.22) 

so that z is indeed singular in this case. L'is = implies that r( m ) and p( m ) 
are parallel; for this reason we call such singularities m-radial. 

L 2 r s = implies that L 2 ^ = for all 2 < j < m, and therefore that 
dJj^-\ = for 2 < j < m. Thus, m-radial singularities have corank at least 
m — 2, and only 3-radial singularities can have corank one. 

In fact, 3-radial singularities are necessarily degenerate. From (J6.22j) . 
we may take C3 = 1 to be the only nonzero coefficient in ()6.8)1 . so that 
K = J - (£,%■.)". Since L? 3 s = r 2 3 ^p 2 3 ^ — (r( 3 ) -p(3)) 2 , K has a single nonvanishing 
6x6 block given by 

( -r (3 )P( 3 )P± r 2 3) P ± 

V -P 2 (3) F ± ^(3)P(3) P ^ 

where is the projection onto the plane in M. 3 perpendicular to r( 3 ) and 
P(3). Since = Pj_, it follows that 

K 2 = 0. (6.24) 

Since 3-radial singularities are degenerate, the singularity formula for the 
Maslov index cannot be applied to them. However, as we now show, they 
do not contribute to the Maslov index. Without loss of generality, we may 
assume that n = 3. Then a 3-radial singularity is of the form 

r = an, p = on, (6.25) 



(6.23) 
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where n is a unit vector in IR 3 . To ensure that the 3-radial singularity has 
corank one, we require that no other singularity conditions are satisfied. That 
is, we assume that a and b do not both vanish, that ()6.9j) is not satisfied 
(ie, (r, p) is not a spherical singularity), and that n is not perpendicular or 
parallel to the e 3 , the unit vector along the 3-axis (ie, (r, p) is neither a 
3-axial singularity nor a 12-radial singularity (see below)). Without loss of 
generality, we assume that a ^ (otherwise, reverse a and b in what follows). 
Let v be a unit vector orthogonal to n, and let w = n x v. Displacements in 
p along v and w with r held fixed are transverse to the 3-radial singularity. 
Therefore, the closed curve C, given by 

r(s) = an, p(s) = on + eu(s), (6.26) 

where 

u(s) = cos27rsv + sin27rsw, (6.27) 

encloses the singularity. For sufficiently small e, no other singularities are 
enclosed. 

To calculate the Maslov index of C, we use the formula (|3.2|) explicitly. 
Letting A, B and C denote the three columns of M, we get that 

. dH m A 

A : = — + 1— = un + e£(s)u(s), 

dLi2 .dL 12 . „ . „ . 

B := — h i—x — = (a - ib)e 3 x n - zee 3 A u(s), 

op or 

dl? dL 2 

C : = — + i— = 2ea{a - ift)u(s) + 2ie 2 ah, (6.28) 
op or 

where a; is a constant and £(s) is complex and of zeroth order in e. (The ex- 
pression for A follows from symmetry considerations; since H is rotationally 
symmetric, its gradients with respect to r and p must be linear combinations 
of its arguments. The expressions for B and C are obtained from straight- 
forward calculations.) Then 

|M|(z(s)) = (C x A)-B = 2ecua(a-ib)((a-ib)u 3 (s)+ien 3 ) + 0(e 3 ), (6.29) 

where u 3 (s) = e 3 • u(s) and n% = e 3 • n. From (J6.29J) . 

arg |M|(z(s)) = const + arg((a — ib)u s (s) + ien 3 ) + 0(e 2 ). (6.30) 

The quantity (a — ib)u$(s) + ien 3 lies on a ray through ien 3 (n 3 7^ by 
assumption), and therefore has zero winding number. The 0(e 2 ) term can 
be neglected. Thus, 

fi(C) = 0. (6.31) 
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6.2.4 c 2 7^ 0, Cj^2 = 0. (l,2)-axial singularities. 

These are singularities of the form dL\2 = 0, which implies, and is implied 

by, 

n = r 2 = pi = P2 = 0. (6.32) 

We call such singularities (l,2)-axial. From (|6.32|) . it is clear that (l,2)-axial 
singularities are also 3-radial singularities, and therefore have corank at least 
two. They do not contribute to the Maslov index. 

6.3 Maslov indices of rotational actions. 

Let 

L (m) = \j L \ m y 3 < m < n. (6.33) 

It is straightforward to verify that the Hamiltonian flow generated by L( m ) is 
a 27r-periodic uniform rotation of r^) and p(j) in their common plane. That 
is, the orbits generated by Lt m ) are of the form 

f(m)( s ) = cossr (m) + sins— ^ (r 2 m) p (m ) - (r(m) ■ p (m ))r (m) ) , 

P(m)(s) = cossp (m) - sins— ^- (p\ m )*{ m ) - (r (m) • P( m ))p( m )) ,(6.34) 

while the components rj and pi with i > m are left unchanged. The L^'s, 
together with L12, which generates 2-7r-periodic uniform rotations in the 12- 
plane, constitute a set of n — 1 action variables associated with rotational 
symmetry. The remaining action variable is obtained from the Hamiltonian 
H (for example, from the reduction of its flow to the radial phase plane). 

It is straightforward to determine the Maslov indices /X12, /i(3), . . . , fi( n ) of 
these rotational actions. Let (r, p) be a regular point. Consider first the angle 
contour through (r, p) generated by the flow of L 12 . This is a 27r-rotation 
of T( 2 ) and p(2) in the 12-plane. The r-orbit can be contracted to a single, 
nonzero point without encountering any co dimension-two singularities (in 
particular, it is readily shown that the contraction can be performed keeping 
r 2 , p 2 and r ■ p fixed, so that no spherical singularities are encountered). The 
p-orbit can then be similarly deformed without encountering singularities. 
It follows that 

/ii 2 = 0. (6.35) 

Consider next the angle contour through (r, p) generated by L( m ). From 
(|6.34|) . the projection of this contour to the (r\,-,pj)-plane for j < m is a pos- 
itively oriented ellipse about the origin. These m ellipses can be contracted 
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in turn to nonzero points in their respective planes. These deformations can 
be performed keeping r 2 , p 2 and r • p fixed, so that no spherical singularities 
are encountered. With each deformation up to and including j — 3, a single 
j-axial singularity is encountered (where rj = pj = 0). Since the ellipses are 
positively oriented and j-axial singularities are elliptic, each contributes +2 
to the Maslov index. The remaining ellipses in the 1- and 2-phase planes can 
be contracted without encountering any singularities. Thus, 

/i (m) = 2(m - 2). (6.36) 

In particular, fj,r n ), the Maslov index associated with the total squared angular 
momentum in n dimensions, is 2{n — 2). This leads via (jl.lj) to the semiclas- 
sical quantisation condition L 2 = (I + (n — 2)/2) 2 , which in turn agrees with 
the exact eigenvalues of the Laplacian on the n-sphere, 1(1 + n — 2) (see, eg, 
jH]) up to an additive /-independent constant. (For n — 3, the semiclassical 
and exact eigenvalues are {I + \) 2 and /(/ + 1) respectively). 

7 Discussion 

For integrable systems in M 2n , the Maslov index of a closed curve in the 
regular component is, under certain genericity conditions, given by a sum of 
contributions ±2 from the corank-one nondegenerate singularities enclosed. 
The sign depends on the stability of the degeneracy and the orientation of the 
curve. We also obtain expressions for the transverse Liapunov exponents of 
corank-one singularities. The fact that the index is unchanged through local 
bifurcations implies relations amongst the stabilities and orientations of the 
singularities involved. For 5'0(n)-invariant systems, we recover the fact that 
the Maslov indices associated with L(j) , the magnitude of angular momentum 
restricted to the first j components, is 2(j — 2). Natural extensions of these 
results would include general cotangent bundles and higher Maslov classes 
|2*H 12*2*] . for which the sources should be (nondegenerate) singularities of 
corank greater than one. 
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